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■ Abstract. High-latitude laminar confinement of the interior field Bi is shown to be possible. Mean 

Q^l downwelling U as weak as 2 x lO^^cms^' - gyroscopically pumped by turbulent stresses in the 

D ■ overlying convection zone and/or tachocline - can hold the field in advective-diffusive balance 

' within a confinement layer of thickness scale 5 ^ 1.5 Mm ^ 0.002/?©. The confinement layer sits 

at the base of the high-latitude tachocline, near the top of the radiative envelope and just above the 
'tachopause' marking the top of the helium settling layer. A family of exact, laminar, frictionless, 
axisymmetric confinement-layer solutions is obtained in cylindrical polar coordinates, for uniform 
downwelling in the limit of strong rotation £2i and stratification A^. The downwelling cannot 
penetrate the helium layer and must therefore feed into an equatorward flow immediately above 
the tachopause. The retrograde Coriolis force on that flow is balanced by a prograde Lorentz force 
Q . within the confinement layer. Buoyancy forces keep the tachopause approximately horizontal. For 

V-j ' typical solar values ^ 10^^ s^' this type of dynamics holds over a substantial range of colatitudes, 

C/3 , e.g. nearly out to colatitude 40° when U < lO^^cm s^' for modest |Bi | values ^ tens of gauss. 

^ ^ ■ The angular-momentum budget implied by the downwelling and equatorward flow, importing 

low and exporting high angular momentum, dictates that the confinement layer must exert a net 
retrograde torque on its surroundings through laminar Maxwell stresses. Some of that torque is 
^ ■ exerted downward through the tachopause upon the interior, against the Ferraro constraint, and 

f*--- I the rest is exerted across the periphery of the confinement layer at some outer colatitude < 40°. 

i The profiles of velocity and magnetic field within the confinement layer are fixed by two external 

' conditions, first the partitioning of the torque between the contributions exerted on the interior and 

across the periphery, and second the vertical profile of Maxwell stress at the periphery. In default of 
i detailed models of what happens near the periphery, we provisionally suggest that a natural simplest 

' choice of model would be one in which all the net torque is exerted on the interior. 
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X : INTRODUCTION 

^ 6 1 

The near-rigid rotation Q.^ = 2.7 x 10 rad s observed in the Sun's interior can be 

most credibly explained via the Ferraro constraint from a confined global-scale interior 

magnetic field Bi (Gough & Mclntyre 1998, hereafter GM98); also Mclntyre (1994- 

2007). For stability must have comparable toroidal and poloidal components (e.g. 

Braithwaite & Spruit 2004). Bj could be axisymmetric and aligned with the Sun's 

rotation axis as proposed in GM98, or oblique as is typical of Ap stars. We focus on 

the aligned case as presenting, in some ways, the greatest problems. 



The main problem, previously addressed by Garaud (2002-08) and by 'B run & Zahn 



( 2006h . is how to confine Bi at the pole and in high latitudes. It is necessary to stop the 



poloidal field from diffusing up through the polar caps and thereby imposing the convec- 
tion zone's high-latitude differential rotation upon the interior. Such differential rotation 
conflicts with the helioseismic evidence. 



FIGURE 1. Streamlines of the downwelling flow u (dashed curves) and of the poloidal part of B (solid 
curves), spacing arbitrary. The horizontal axis is colatitude in arbitrary units, measured from the central 
symmetry axis or rotation axis. The vertical axis is heliopotential altitude, or radial increment, z, in units 
of 5. Although there is a field line on the central symmetry axis, as well as a streamline, the strength of 
the field decays like exp(— with altitude z- The downwelling profile m,(z) — —U[l — exp(— 3z/5)]. 



Following GM98 we propose that Bj is confined to the interior, in high latitudes, via 
a frictionless laminar boundary layer at the base of the stably stratified tachocline. The 
confinement is brought about by weak downwelling of magnitude U , say, taken to be a 
persistent feature of the mean meridional circulations (MMCs) in the lower tachocline. 
Such downwelling is to be expected from the 'gyroscopic pumping' cause d by turbulent 
Reynolds and Maxwell stresses in the overlying layers (|McIntyre l2007h . in the same 
way that Ekman pumping is caused by ordinary viscous stresses. 

The ove rlying turbule nt layers consist of the convection zone and possibly also the 
tachocline (|Spruitil2002r) . There are uncertainties in how to characterize those turbulent 
layers in detail. But because of the 'antifrictional' sense of the convection zone's tur- 
bulent stresses - driving it retrogradely relative to the interior in high latitudes - there 
should be a systematic tendency, by one means or another, for the gyroscopic pumping 
to produce downwelling above the to p of the radiative envelope i n high l atitud es. One 
possible such scenario is discussed in lMcIntyre (l2007l) . following [SpruS (l2002h . There 
the gyroscopic pumping is, in fact, a case of ordinary Ekman pumping near the bottom 
of an MHD-turbulent tachocline. 

We assume that the bulk of the radiative envelope beneath is itself locked into rigid 
rotation with the help of the Ferraro constraint from Bi, and that gravitational set- 
tling has produced molecular- weight gradients in the form of a helium settling layer in 
the outer 50-1 0() Mm o r so. Such a helium layer is a feature of standard solar models 
( e.g. Ciacio et al.l 1997 ). Its existence is indicated also by the helioseismic evidence, de- 
spite current uncertainties about heavy-elemen t abundances and their effects on opacity 
(|e.g. Christensen-Dalsgaard & Pi Maurdl2007|) . 

The helium diffusivity is tiny, ~ 10cm'^s~\ in comparison with the thermal and mag- 
netic diffusivities K ~ lO^cm^s^^ and T] ~ 3 x l O^cm^s"^ (e.g. Gough 2007). Therefore 
the helium layer is nearly impervious t o MMC s. Mestel & MossI (|l986l) call this the '/i- 
choke' or '/i -barrier' effect; see also ^ Mestel (11953'). The high-latitude downwelling, 
whatever its origin, must therefore feed into an equatorward flow just above the 'ta- 
chopause' marking the top of the helium layer. The retrograde Coriolis force on that 
equatorward flow needs to be balanced by a prograde Lorentz force. Buoyancy forces 
from the stable stratification keep the tachopause and the helium layer beneath it very 
close to the horizontal, along with the stratification surfaces themselves (Mclntyre 2007, 
§8.5). We report a new family of exact steady solutions of the nonlinear equations show- 
ing how all these elements fit together, confining B; within a layer of thickness scale 
d = rj/U while transmitting a retrograde torque to the interior. 




FIGURE 2. Further views of the same solution. The view (a) of field lines shows only those lines that 
emerge from the tachopause at a single colatitude near the pole. The lines emerging at other colatitudes 
have the same shape dilated horizontally (see text). The profiles (b) are those of (solid), Ur (dashed) 
and (dot-dashed), in units of U, Ur/ 5 and 5QUr/ 5 respectively; (c) are (solid), Br (dashed) Bij, (dot- 
dashed), and the auxiliary function B,j, (dotted), in units of (rjil;)'/^ for B- and (T7£2i)^''^r/5 for the rest. 
The shape of (z), not shown in (b), is nearly indistinguishable from u,j, (z) — u,j, {°o). See dUi-® below. 

CONFINEMENT-LAYER SOLUTIONS 

The confinement-layer solutions are axi symmetric similarity solutions obtained in cylin- 
drical polar coordinates (r, 0, z) and valid in some region surrounding the pole, for any 
finite stratification A'^. Figure [H shows an example, with the north pole at the centre. The 
z axis is central, pointing upward, and z measures altitude in units of 5. The cylindrical 
radial coordinate r is proportional to colatitude and is in arbitrary units, the structure be- 
ing self-similar under radial dilatation. The similarity solutions have vertical field com- 
ponents Uy and By that are independent of r, and horizontal components Ur, and Br, 
Bp proportional to r. Because V ■ B = we have 2Br = —Bir where the prime denotes 
d/dz. We assume anelastic flow with 8 <^ pressure scale height (60 Mm or more), so that 
background density is constant (Boussinesq limit) and V ■ u = 0, implying 2ur = —uir. 

The boundary conditions include continuity of B at the tachopause z = 0. The mag- 
netic field Bi just beneath the tachopause therefore has a matching similarity structure. 
The simplest such field has Bir — Bior and B^y = —2B[qz + Q, where B^q and Q are con- 
stants, so that V ■ Bi = 0, while the azimuthal component 5;^ is taken oc Bi^r to ensure 
the vanishing of the interior's azimuthal Lorentz force, right-hand side of ([3]) hereunder. 
Figure [2l shows further views of the same solution. The equations solved are 

riB'^ = B[u,-B,u[ (1) 
r\B'^ = B'^Uz-B^u'^ (2) 

2a,ur = -a^u'^r = B,B'p-B'fy = Bl{Bp/B,)' (3) 

of which the first two come from the induction equation after substituting the similarity 
structure into V x (u x B), and the third from the azimuthal momentum equation in the 
limit of small Rossby number, thus retaining only Coriolis and not relative accelerations; 
B is in units of Alfven speed, ~ ^cms^^ per gauss at tachocline mass densities. The T] 
terms are exact because the similarity structure makes the horizontal derivatives vanish. 



In the absence of specific information about the vertical profile of Maxwell stresses 
BrB(j) at the periphery, there is an undetermined function of z in the problem. We may use 
this freedom to specify the vertical profile u^iz) of the downwelling. Then ([T]) becomes a 
linear problem for B^{z) ■ It is solvable with B'^{0) = — 25io (matching the r components) 
and with B^{z) decaying like exp(— as z ^ 0° . This determines both ^^(z) and Q. 
Then, provided only that |%| converges to U faster than exp(— 2z/5), © has a solution 

u 

^dz. (4) 

In the example shown in the figures we took = — t/[l — exp(— 3z/5)], which implies 
that Ur = — ju'^r = l(Ur/5)exp{—3z/5), giving the 5^ profile shown as the dotted, 
leftmost curve in Fig.l2t. Then (O has a corresponding solution 

f^'nB''-u,B' 
= u<p = / dz. (5) 

Jz B^ 

Under our assumptions, both @) and dS]) are finite at z = 0, and are evanescent as z ^ 00, 
respectively like exp(— 2z/5) and exp(— z/5) in the example shown. The foregoing 
procedure is robust and well-conditioned. In the simplest cases in which and u'^ are 
both negative for all z, as in the figures, it is easy to see by inspection of qualitative- 
ly integrating it downward from z = °°, that 5" on the left and both terms on the right 
are positive for all z. It then follows from dH) that Btp is negative for all z. But m,^ and 
can change sign, though in fact a,^ is negative and monotonic in the example shown. 

We still have a pair of undetermined parameters at our disposal because ^ and Q 
are also satisfied, for any constants a and /3, by 

B(j, = Btp + aB^r and Utp = + (ocmj- + /3)r . (6) 

Here ([U) has been used to simplify the last term on the right. The a term in B^j, contributes 
nothing to the azimuthal Lorentz force on the right-hand side of (O, but does change 
the Maxwell stresses acting across the tachopause z = and the periphery, r = say, 
by equal and opposite amounts. In other words a governs the partitioning of Maxwell 
torques between tachopause and periphery. 

In the example shown a was chosen, purely on Occam's-razor grounds, to make 
the Maxwell torque on the periphery zero. Lacking information about conditions at the 
periphery, zero is arguably the simplest choice. The torque <x ByB^ dz. It is zero if 

a=-f||4 = 1.315-' (7) 

Jo BrB^rdz 

from the first of ^\ the quotient is independent of r because of the similarity structure. 
To find /3 we integrate © across the tachopause and use Q to give 

u^{0) = -A-^Ur{0) (8) 

where A = Bj{0) / {2Q.iri), the Elsasser number based on By{0) (~ BioS), determining 
the direction of the frictionless slip flow at the tachopause just above the rigidly -rotating 
interior. The slip flow is equatorward and retrograde, following a logarithmic spiral. 
For finite viscosity v there would be a laminar Ekman layer of thickness scale 5v ~ 



y/^y/^ ~ 3 X 10"^ Mm < 5, if we take v ~ SOcm^s^i dOoughl llOOTh . Its flow 



is unobstructed by the field lines since magnetic diffusion on the scale 5v is almost 
instantaneous. Note ([8]) gives the spiral just above the Ekman layer. From ^ and ([8]), 

/3 = -r-^u^{Q) - A-^r-^Ur{Q) = 1.73t/5"^ x 10^ . (9) 

In fact A^^ measures the spiralling of the field lines as well as that of the flow lines, 
because dH)-® imply the order-of-magnitude relations 

5^ ~ 5^ ~ laiUr/B, , ~ ~ UB^/B, , (10) 

and 5^/5, ~ m^/m, ~ . (11) 

Recall that 5 = rj/U and 2ur = —u[r ~ Ur/5, 2Br = —Bir ~ Bior. The numerical 
factors implicit in (fTTl) differ considerably from unity because of the disparity in verti- 
cal scales between exp(— z/5), exp(— 2z/5) and exp(— 3z/5), along with the peculiar 
balance of terms in © that enables B^j, to evanesce faster than exp(— z/5). In the ex- 
ample shown in the figures, A = 1.57 x 10^^. Fig. |2]shows that attains values 
considerably smaller numerically, and \U(j)/ur\ distinctly larger, than A^^ ^ 60. 

By contrast with GM98's thermomagnetic boundary layer, whose dynamics crucially 
involved the tilting of stratification surfaces, our exact solutions of ([I])-® impose no 
restriction on U values and mass throughput for given B^q. However, there is an implicit 
restriction, for given peripheral radius r = rp and stratification A^. has been assumed 
strong enough to hold stratification surfaces horizontal. Only then can the uniform 
downwelling satisfy the thermal diffusion equation, displacing the stratification surfaces 
vertically without tilting them. A scale analysis, omitted for brevity, shows that the tilting 
becomes noticeable at the periphery if in order of magnitude 

(a5io/rp)i/3, b{B,orp)-' 



U ~ C/crit = min 



(12) 



where a = OA{ri /Ky/^{ri^N/Qi) ~ 0.7 x lO^cm^s^^ and b = 0.1{ri /xy/^riN ~ 

0. 15 X lO^cm^s^. The min function arises from the azimuthal vorticity balance. In 
the strong-field case (second argument, roughly corresponding to A > 1), the tilting of 
the stratification surfaces is balanced solely by a Lorentz force-curl. In the weak-field 
case (first argument, A < 1) there is an additional contribution from vortex twisting 
2Q.iu'^. The crossover corresponds to B^orp ~ 15 cms^^ (~ 30G) when rp = 350Mm, 

1. e. to |Bi| ~ 30G near an outer colatitude ~ 40°. Then Ucm ~ lO^^cms^^. 

The Rossby number M = maxlr^^u- V(ra^)/(2aiM^)| ~ {U/B^{0))^ ~ Uy{riBio)^, 
< 10-6 at crossover ii U < Ucm, and similarly small throughout the parameter 
range of interest, strongly justifying our use of the smaW-M limit. Inverse gradient 
Richardson numbers \u'\^/N'^ ~ (U /Ucni){r} /k)M in the confinement layer, and 
~ {U/Ucrit){r]/K)M{5/5y)^ ~ {U/Ucnt){vMU^^i/{vBl) in the Ekman slip layer, 
re spectiv e ly 5, 1 0"^^ and < 10^^ at crossover. So the flows are strongly shear-stable. 



Spruiti (| 19991) shows that in stably stratified shear flows the first MHD instabilities 
to kick in will be diffusion-mediated Tayler kink or tipping instabilities of 5^. In the 
weak-field case A < 1 (the most vulnerable, with strong spiralling) we find Tayler 
stability for i7 < (5io/'BiOcrossover)^'^^f^crit- Stability increases further when A > 1 . So 



the solutions probably represent real laminar flows. 

CONCLUDING REMARKS 

The suggestion in (fT2l) of a limiting mass flow and therefore, by implication, of an 
upper bound on the torque transmissible to the interior, is no more than a suggestion 
at present. However, the scaling leading to (fT2l) does have points of similarity to the 
scaling governing the mass-flow-limited thermomagnetic boundary layer proposed in 
GM98. A peripheral thermomagnetic boundary layer might impose a mass-flow limit. 
Such a limit would have implications, in turn, for the possible range of interior field 
strengths |Bi|. In particular, the steep falloff of Ucm on the strong-field side of (fT2l) 
suggests a sharp upper bound on confinable |Bi| strengths. 

A mass-flow limit, if confirmed, would also bear on the question of whether a Tayler- 
Spruit dynamo can run continuously or intermittently in the tachocline above the con- 
finement layer (Mclntyre 2007, §8.4). That question is critical to associated questions 
about deep tachocline ventilation and lithium burning. 
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